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In the framework of the minimal super symmetric standard model (MSSM) and 
the R-parity breaking model (RBM), we investigate production processes of the scalar 
top (stop) at HERA energies. These models are characterized by the possible exis- 
tence of the light stop whose mass is lighter than those of the top quark and the other 
squarks. It is shown that in the MSSM the stop pair production via boson-gluon fu- 
sion gives a sizeable cross section and the most serious background ep — > eccX could 
be suppressed by appropriate kinematical cuts. We also show that in the RBM the 
stop is singly produced in the neutral current processes and we have a clear signal as 
a sharp peak in the Bjorken parameter x distribution of the scattered electron. 



I. INTRODUCTION 

The standard model (SM) of strong and electroweak interactions describes amazingly well 
particle physics phenomena up to the electroweak breaking scale. Particularly, the agreement 
between theory and experiment in the LEP experiments was so impressive. Nevertheless, 
we know that the SM could not be the ultimate theory of everything, so long as there are 
a large number of free parameters, the arbitrariness of particle masses and mixing angles, 
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and the lack of any explanation for the replication of generations and so on. Still there is 
no experimental sign of new physics beyond the SM. 

Various theoretical attempts are open to go beyond the SM. Among them the super- 
symmetry (SUSY) |IJ seems to provide us an elegant and reasonable approach to meet the 
expectation, provided that the SUSY breaking scale is O(lTeV). In the framework of SUSY, 
every boson (fermion) is accompanied by a fermionic (bosonic) partner with equal mass and 
coupling. The requirement of equal mass is relaxed by the broken SUSY either spontaneously 
or explicitly through introducing soft SUSY breaking terms like gaugino mass terms, scalar 
mass terns or trilinear couplings of scalar particles. Still, there is no experimental evidence 
for SUSY partners of any of the known particles, and only some limits on their masses are 
given . The theoretical expectation for SUSY, however, is getting stronger than ever. The 
analysis of the running coupling constant by SUSY GUTs recently done seems to encourage 
people who believe in SUSY [||]. 

The purpose of the present paper is to focus our attention on the scalar top quark (stop 
ti) and investigate its production mechanisms and decay properties in the framework of the 
minimal supersymmetric model (MSSM) and the R-parity breaking model (RBM) [Q. We 
calculate cross sections together with detectable signals expected to observe at HERA, the 
world first electron-proton collider. Contrary to predicted high mass of the top quark, which 
is not discovered yet, the stop mass could be lower than that of the top quark in a model ||. 
The expected mass of the stop is clearly within the reach of HERA. This is one of our strong 
motivation of the present work. Higher energies available than LEP and cleaner events than 
hadron colliders are another reason which has driven us to HERA. Particularly, the stop 
could be singly produced in ep collisions in the RBM ||[7],||] because of the existence of the 
electron-quark-squark couplings. Clearly HERA is the best machine to search for the stop 
in such models, because we expect the remarkable peak structure in the Bjorken parameter 
x distribution. The theoretical prediction of the discovery limit of the stop depends upon 
various model parameters as well as the energy and luminosity of HERA. We search for a 
possible experimental signal observable in ep collisions. 
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The present paper is organized as follows. After the model building of the stop in 
Sec. II we discuss the stop production in MSSM in Sec. III. Section IV is devoted to the 
resonance production of the stop by invoking the R-parity breaking. Finally, in Sec.V we 
give concluding remarks. 



II. MODEL BUILDING 



In the framework of the supersymmetric standard model, scalar fermion mass matrices 
are expressed by 

/ \ 

m~ atmt 

M}= '* , (1) 

V ^ m 7n J 

where m j L R an d «/ are the SUSY mass parameters and my denote the ordinary fermion 
masses. We can see from Eq. (1) that for the sleptons and the squarks except for the stops, 
the left and right handed sfermions are mass eigenstates in good approximation owing to 
small fermion masses in the off-diagonal elements of the mass matrices. On the other hand, 
the large mixing between the left and right handed stops will be expected from the large 
top-quark mass ||, and the mass eigenstates are expressed by 
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We can easily calculate the mass eigenvalues of the stops as 
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The SUSY mass parameters mj L and at are parametrized in the following way 
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= M| + m| cos 2/3 ( - - - sin 2 W ) + m 2 , (5) 



/ , M 4 + g m i cos 2 ^ sin °w + mi, (6) 
a t = -(A t + fi cot (3), (7) 

where Mj , tan /3, /x and A t denote the soft-breaking scalar masses, the ratio of two Higgs 
vacuum expectation values (= V2/V1), the Higgs mass parameter and the trilinear coupling 
constant, respectively. In Fig. 1 we show contours of the various lighter stop masses 
in (A t , M~) plane for given values of (mt, A 4 ; tan/5), where we take Mj l =Mj r =Mj. We 
find that if SUSY mass parameters and the top mass are the same order of magnitude, the 
cancellation could occur in the expression for the lighter stop mass Eq. (f|). So we get one 
light stop t\ lighter than the top as well as the other squarks for a wide range of the SUSY 
parameters. 

After the mass diagonalization we can obtain the interaction Lagrangian in terms of the 
mass eigenstate t±. In Fig. 2 we show the Feynman rules for the lighter stop t\ couplings 
to the gauge bosons 7, Z and g. While the stop couplings to the gluon and photon do 
not depend on the mixing angle 8 t , the Z-boson coupling depends sensitively on 6 t . More 
specifically, it is proportional to Cj = | sin 2 9\y — \ cos 2 6f Note that for a special value of 
6' t ~0.983, the Z-boson coupling completely vanishes . 



III. STOP PRODUCTION IN MSSM 



In the MSSM the light stop t± will be produced in pair via boson-gluon fusion ||11|| : 

e ± p -> e^JlX. (8) 
Feynman diagrams for the sub-processes are depicted in Fig. 3. Referring to Fig. 3 we find 



that the 7 — g contribution to the sub-processes is much larger than that of Z — g |TT 
Since the 7 — g contribution does not depend on the stop mixing angle &t, the total cross 
section is not sensitive to 9t, but only depends on the stop mass mr . Owing to the photon 
dominance, we can use the Weizsacker- Williams approximation (WWA) ; 



^ = ldyP(y)JdrjG(rj,s')^ (9) 

where y — {p ■ q)/{p ■ £ e ) and z = {p ■ Pf)/(p • q) with the variables defined in Fig. 3. The 
WWA factorizes the cross section of the process (|8|) into the probability for emitting photon 
from lepton ; 

P(,) = iil±(l^)!,o g «^, (10) 
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the boson-gluon fusion cross section involving the real photon 
da 4 naa s 



2mf - 2mfs'z(l - z) + s V(l - zf (11) 



dz 9J'¥(l-2) 2 

and the gluon distribution function G(r],s') \T2\, where rj denotes the momentum fraction 



of the gluon in the proton and s' = yr]s. The numerical integration has been performed 
by using the program package BASES We have the total cross sections : 2.4pb for 

m T = 30GeV, 0.5pb for m T = 40GeV, 5 x 10^ 2 pb for m T = 60GeV and 8 x l(T 3 pb for 
m T = 70GeV. 

Now, what is the experimental signature ? How to suppress its background ? Actually, 
the stop can decay into the various final states : 

t 1 ^tZ 1 (12a) 
^bWi (12b) 
^btv (12c) 
^bvt (12d) 
^bWZ x (12e) 
^bivZ x (12f) 

(12g) 

where Zi, W\, v and t, respectively, denote the lightest neutralino, the lighter chargino, 
the sneutrino and the slepton. If we consider the light stop with mass lighter than 40GeV, 
the first five decay modes (12a) to (12e) are kinematically forbidden due to the lower mass 



bounds for respective particles ; m t ~90GeV, m% ~20GeV, m^~45GeV, m^~45GeV and 
m~~40GeV 0. So there left (12f) and (12g). The one-loop mode t\ — > cZ\ (12g) domi- 
nates over the four-body mode t\ — > bivZ\ (12f ), because of the large logarithmic factor 
\n(M p i anc k / mw) appearing in (12g) [H. So we can conclude that such light stop will decay 
into the charm quark jet plus the missing momentum taken away by the neutralino with 
almost 100% branching ratio. 

Naively, it will be expected that TEVATRON and/or LEP can set severe bounds on 
the stop mass through the processes ; gg —+ t\t\ — > ccZ\Zi (TEVATRON) and/or Z — > 
till (LEP). However, the situation is not so obvious. Baer et al. Jl4| have performed the 
analyses of the experimental data of 4pb _1 integrated luminosity TEVATRON running, and 
have obtained the results that the stop could easily be escaped the detection if ~ 
lOGeV. Such large neutralino mass could make the charm quark jets softer. Consequently 
the stop production cross section plotted against the missing transverse energy becomes 
smaller than the present upper bounds, where we impose cuts on the missing transverse 
energy. Moreover, we find that LEP cannot exclude the light stop for appropriate mixing 
angle Qf In Fig. 4 we show the stop mass dependence of the decay width T(Z — > tit*). The 



horizontal line corresponds to the present upper bound for ATz at LEP ||15|| . We can find 
that the stop with mass mj ~ 20GeV has not been excluded if the mixing angle 6 t is larger 
than about 0.6. The origin of such sensitivity of T[Z — > titj) has already been discussed 
in the previous section. Recently DELPHI has performed the analyses on the direct search 
for the stop through e + e~ — > t\t\ and also studied ATz |15| . The result is still negative. 
According to their analyses we have the excluded region in (8 t , m^j) plane shown in Fig. 5. 
In the same figure we show the excluded region at TRISTAN assuming massless photino 
fl6|j . In this figure we also show the typical searchable regions at HERA. We can clearly see 
that HERA would cover the window near 6 t ~1 in LEP region. The upper bound at HERA 
goes upward for large integrated luminosity and efficient background suppression. 

For the light stop (m^~40GeV) production at HERA, the event signature is eccX with 
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large missing momentum due to unseen neutralinos. The most serious background process 
comes from the direct charm pair production ep — > eccX via boson-gluon fusion process 
whose total cross section is about 1/zb without any cuts. In order to suppress this background 
we find the following kinematical cuts to be efficient ; l)the lower transverse momentum 
P^-^ cuts for both D and D mesons in laboratory (LAB) system, 2)the upper momentum 
transfer squared Q 2 cut and 3)the lower missing transverse momentum^ cut in LAB system 
4)the acoplanarity angle 4>(d-d) cu ^' 4>(d-'d) Dem g the angle between the D-proton and D- 
proton planes in the boson-proton center of mass system. 5)the polar angle d -q cuts, D ^ 
being the polar angle of D or D with respect to the initial proton axis in LAB system. 
In Fig. 6 we show the acoplanarity angle distribution for the stop production followed by 
the decay t\— >cZ\-^ (D-meson) plus missing momentum and those for the background with 
cuts mentioned above. These events in Fig 6 were generated in the kinematical regions ; 
5<Q 2 <100(GeV/c) 2 , 10~ 3 <x<l, 0.1 < \sin9 D ^\ andP T > 5GeV/c by making use of the 
generator AROMA with JETSET |Tj§. Here we took = 30GeV, m ? = 18GeV and 
the integrated liminosity L = 50pb _1 . We see that the D-meson jets with cos^^ ^^O 
can be used to separate the stop production from its background coming from boson-gluon 
fusion process. 



IV. STOP RESONANCES IN RBM 

In the RBM there could exist the exotic Yukawa couplings of the stop ti, which violate 
the R-parity defined by (_) 3B + L + 2S . Here L, B and S denote the lepton number, baryon 
number and spin, respectively. The R-parity breaking (RB) couplings in the supersymmetric 
models |JJ] are required in order to explain the cosmic baryon number violation, the origin 
of the masses and the magnetic moments of neutrinos and some interesting rare processes 
in terms of the lepton and/or baryon number violation. Here we take a coupling of the stop 
h ; 

Aat = A' 131 COS OtihdPLB + t\eP R d) (13) 
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originated from the RB superpotential W = X'^LiQjDk, where i ~ k are generation indeces 
and L, Q and D donote the chiral superfields. The couplings generated from the super- 
potential are most suitable for the ep collider experiments at HERA because the squarks 
will be produced in the s-channel in e-q sub-processes. Here we consider the coupling ([T5D , 
which generates the s-channel stop t\ production in the neutral current (NC) process ||: 

e ± p (hX) -> e ± qX (14) 

and the relevant Feynman diagrams are depicted in Fig. 7. The production of squarks in 
the first and second generation at HERA has been discussed in Ref. 

We calculate the inclusive differential cross section for the NC processes e ± p — > e qX 
with polarized beams : 

^Q2^U = 5prE [^Q 2 )T,Ti{ei R q) + g(s,Q a )£T,( e £ B ff)] , (15) 

where x and Q 2 denote the Bjorken scaling parameter and the four momentum transfer 
squared, respectively, and q(x, Q 2 ) are the quark distribution functions ]T2[] in the nucleon. 



The analytic expression for the coefficients Tj(e^ R q) is presented in the Appendix. The cross 
section depends sensitively on the decay width of the stop. In this calculation we assume 
BR(ti -> ed)~100%, i.e., 



a 



r~ = -FRBCtOr^, (16) 
where the coupling strength -Frb^i) is defined as 

ffcefc) = (17) 

This situation corresponds to the case of <m t + , m& + . While for >m t + 
or m^ + m^ BR(£x — > ed) will compete with BR(ti — > tZi) or BR(tx — > bW\). The analyses 
for such case will be presented elsewhere [|T!)J . 

We expect a clear signal of the stop as a sharp peak in the Bjorken parameter x distri- 



bution. The peak point corresponds to x — m~ / s. In Fig. 8 we show the x distributions 
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at fixed Q 2 . As clearly seen from Fig. 8, the lower Q 2 cuts would be very efficient for 
suppression of the SM background, since the s-channel stop contribution is independent of 

Q 2 . 

In Fig. 8 we take the beam to be e~. The event rates depend on not only the RB 
coupling strength A' 131 but also the kind of beam. This is shown in Fig. 9, where we fixed 
x = 0.2 at the peak value. It is found that the e + beam is more efficient than the e~ one to 
separate the stop signal from the SM background. This can easily be understood from the 
structure of the coupling. While the e~ collides only with sea <i-quarks in the proton, the 
e + collides with valence ci-quarks. The difference of the structure function of the proton is 
naturally reflected in the cross sections . We also note that the polarized and beams 
could be advantageous to suppress the SM background. In Fig. 10 we show the y{= Q 2 / sx) 
dependence of the asymmetries defined as 

c = da(e+)/dxdy - daje^/dxdy 
dcr(e^) /dxdy + da{e^)/dxdy 

and 

A _ = da(el)/dxdy - daje^/dxdy 
da (e^)/ dxdy + da (e^)/ dxdy 

We find that the polarized e ± beams will be effective to identify the stop signal. Next we 
show in Fig. 11 the searchable parameter region at HERA in (A' 131 , mqr) plane. The shaded- 
region is experimentally excluded through the atomic parity violation experiments [|]]. The 
area inside the solid contour represents the region accessible at HERA for the production of 
more than ten signal events above the SM background with Q 2 >10 3 (GeV/c) 2 in 100pb _1 
running. HERA will open up a large window as for the existence of the stop with RB 
coupling at mass up to 200 (270) GeV at the coupling of A' 131 ~ 0.1 with e~ (e + ) beams. 
It is well known that the similar x peak could be expected in the leptoquark production at 



HERA PD" |. We should point out that the stop with the RB couplings could be discriminated 
from most of the leptoquarks by its distinctive properties ; l)the x peak originated from 
the stop would exist only in the NC (not exist in the CC) process because there is no RB 



stop couplings to the neutrinos, 2)the e + beams are more favorable than the e~ beams as 
mentioned above. One of the leptoquarks S%/2 with the charge Q = —2/3 will give same 
signature with the RB stop, if the stop has BR(ti — > e<f)~100%. Again we note that this 
situation corresponds to the case of wq x <mt + , nib + m w 1 - While for m^>m t + m% 
or m b + rrifir BR(t! — > tZ x ) or BRfTx — > bW\) will compete with BRf^ — > ed) and the stop 
could be discriminated from the leptoquark Si/ 2 with the charge Q = —2/3 



V. CONCLUDING REMARKS 

We have investigated various production processes of the stop at HERA energies in 
the framework of the MSSM and the RBM. For the stop with mq ~50GeV appropriate 
kinematical cuts will enable us to see signals well above the background. The stop produced 
via R-breaking interactions shows a sharp peak in the x distribution due to its s-channel 
resonance. 

The discovery of the stop would reveal us simultaneously the existence of the top flavor 
and the supersymmetry. This is a really big impact to the present day particle physics. 
HERA would be the ideal machine for our stop searches because it is ep collider with just 
suitable energy range for our purpose. 

Thanks are due to a Grant-in- Aid for Scientific Research of the Ministry of Education, 
Science and Culture(tj03240277). 

APPENDIX A 

The analytic expression for the cross section Eq.(|T5|) is given as follows. 

£ If [el ii? ] Q 2 )u(x, Q 2 ) + £ f^elj (x, Q 2 )u(x, Q 2 ) (Al) 
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for e beams. The formula for e + beams can be obtained by the following replacement in 
the above formula for the e~ beams ; 



q^q 



(A21) 
(A22) 
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(A23) 



Here, denote the electromagnetic charge of matter fermion /, and 



Al = - 1 * (A24) 
cos Ow sin Ow 

A f R = e f tan6 w , (A25) 
where T/ are the third component of isospin and 9w is the Weinberg angle. 
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FIGURES 

FIG. 1. Contour of lighter stop mass m~ (GeV) in (A t , M~) plane. We take M~ =Mj =Mj 
and (mt(GeV), /i(GeV), tan/3)= (150, -300, 2) for (a) and (150, -300, 12) for (b). 

FIG. 2. Feynman rules for the lighter stop t\ couplings to g, 7 and Z. 

FIG. 3. Feynman diagrams for sub-process e ± g — > e ± iii*. 

FIG. 4. m~ dependence of decay width T(Z —> t\t\). The horizontal line corresponds to the 
present upper bound for /STz at LEP. 

FIG. 5. Excluded parameter region in (9 t , mj ) by LEP and searchable region at HERA. Region 
under the dotted line corresponds to excluded region at TRISTAN assuming massless photino, 

which is valid for m~ — m~ > 8GeV. 

ti 7 

FIG. 6. Acoplanarity distribution. Kinematical cuts are > 2GeV/c and 5GeV/c, 5 < Q 2 

< 100(GeV/c) 2 ,^r > 5GeV/c and | sin 6> ^| > 0.1. 

FIG. 7. Feynman diagrams for sub-process e ± q — > e ± q. 

FIG. 8. x distribution at fixed Q 2 . Adopted parameters are m~ = 200GeV, A' 131 = 0.25 and 
9 t = 0.0. 

FIG. 9. y distribution at fixed x using the electron and the positron, x is fixed at 0.2. Adopted 
parameters are = 140GeV and Qt = 7r/4. 

FIG. 10. y distribution at fixed x of differential asymmetries Cr(&) and A e -(b). x is fixed at 
0.2. Adopted parameters are mj- = 140GeV and 9 t = 0. 

FIG. 11. Searchable parameter region at HERA in (A' 131 , fnj)- Kinematical cut is 
Q 2 > 10 3 (GeV/c) 2 . 
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